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We  measure  the  plume  temperature  and  species  concentration  profiles  of  a  heated  mixture 
of  ammonia  and  steam  flowing  through  a  miniature  converging-diverging  nozzle  (area  ratio 
5:1,  exit  diameter  4.57  mm).  The  nozzle  operating  conditions  are  intended  to  simulate  the 
plume  of  a  thruster  system  operating  in  vacuum  environments.  Wavelength  modulation 
spectroscopy  (WMS)  is  a  laser  absorption  spectroscopy  technique  that  allows  for  quantitative, 
time-resolved,  sensitive  diagnostics  of  gaseous  flows.  The  goal  of  the  measurements  is  to  assess 
the  potential  of  WMS  as  a  diagnostic  technique  for  characterizing  this  type  of  system.  The 
measurements  use  three  multiplexed  fiber-coupled  tunable  diode  lasers:  two  for  measuring 
temperature  and  H2O  concentration  and  one  for  measuring  NH3  concentration.  The 
multiplexed  laser  beam  is  spatially  translated  across  the  plume  on  a  two-axis  stage.  Using 
Abel  transforms,  we  are  able  to  successfully  resolve  the  structure  of  the  under-expanded 
plume  despite  absorption  signals  as  low  as  9e-4.  As  a  result,  WMS  demonstrates  excellent 
promise  as  a  sensitive  diagnostic  tool  for  resolving  the  profile  of  such  plumes. 
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Nomenclature 

path  averaged  absorbance  area  of  feature  at  wavelength,  X 

mole  fraction  of  species,  i 

incident  intensity  on  gas  sample 

transmitted  intensity  out  of  gas  sample 

pathlength  probed  by  laser  sensor 

wavelength  of  light 

optical  frequency  of  light 

lineshape  function  of  absorption  feature 

linestrength  of  absorption  feature 

temperature  of  gas  sample 

outer  radius  of  jet 

radial  position  within  jet 

vertical  position  within  jet 

X  component  of  lockin  amplifier  at  2/ 

Y  component  of  lockin  amplifier  at  2/ 

root-sum-square  magnitude  of  the  X  and  Y  components  of  the  lock-in  outputs  at  If 
radial  distribution  of  absorbance  area  of  feature  at  wavelength,  X 
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I.  Introduction 

Plume  diagnostics  are  important  tools  for  characterizing  the  operation  of  rockets  and  other  propulsion 
systems.  However,  many  systems  are  difficult  to  probe  due  to  their  size  (both  large  and  small),  high  velocities,  caustic 
products,  and  varying  pressure  conditions.  Laser  based  diagnostics  are  ideal  for  these  environments  as  they  are  time- 
resolved,  non-intrusive,  and  can  be  very  sensitive.  Time -resolved  measurements  of  species  concentration  and  plume 
temperature  can  be  used  to  validate  models  and  quantify  the  state-of-health  of  these  systems.  Wavelength  modulation 
spectroscopy  (WMS)  is  a  method  of  absorption  spectroscopy  particularly  well  suited  for  these  measurements  because 
it  is  quantitative,  reduces  noise,  and  has  improved  sensitivity  over  traditional  tunable  diode  laser  methods  [l]-[2]. 

In  absorption  spectroscopy,  light  is  passed  through  a  gas  sample  and  will  be  absorbed  if  the  wavelength  of  the  light 
is  resonant  with  a  molecular  transition.  For  homogenous  mixtures,  absorbance  (a)  of  the  light  at  a  given  wavelength 
(A),  is  quantified  by  Beer’s  Law  (1).  Io  represents  the  initial  intensity  of  the  light,  and  It  represents  the  transmitted 
intensity  of  the  light  after  it  has  passed  through  the  sample.  From  the  absorbance,  gas  properties  such  as  species 
concentration  (Xi),  pressure  (P),  pathlength  (L)  and  temperature  (through  the  linestrength  SA(T))  can  be  calculated,  if 
the  lineshape  function  ®  is  known  (or  if  an  entire  isolated  lineshape  is  measured,  since  0(A)dA  =  1).  Transition 
specific  parameters  such  as  5^(7)  and  O  are  calculated  from  spectroscopic  parameters  found  in  databases  such  as 
HITRAN  2012  [3]  or  from  experimental  studies  [4]— [6]. 

a{X)  =  Sx(T)PXtL<KX)  =  -In  (j^j  (1) 

Tunable  diode  lasers  (TDLs)  are  commonly  used  in  absorption  spectroscopy  due  to  their  small  size,  low  cost  relative 
to  other  laser  sources,  and  robust  nature  (resulting  from  of  their  development  by  the  telecom  industry).  The  wavelength 
of  a  TDL  tunes  approximately  linearly  with  the  tuning  of  the  laser  injection  current,  allowing  the  center  wavelength 
of  the  TDL  to  be  swept  over  an  absorption  feature  [7].  WMS  is  a  modification  of  tunable  diode  laser  absorption 
spectroscopy  where  the  wavelength  of  the  TDL  is  rapidly  modulated  at  a  frequency  f  (here,  approximately  100  kHz) 
[2] .  In  the  presence  of  absorption,  this  modulation  gives  rises  to  harmonics  of  the  original  modulation  frequency  (f)  in 
the  detected  intensity  signal.  The  second  harmonic  (2f)  is  of  particular  interest,  as  the  magnitude  is  linearly  related  to 
absorbance  and  can  therefore  be  used  to  infer  thermodynamic  properties  of  the  gas.  If  the  second  harmonic  signal 
alone  is  used  for  thermodynamic  measurements,  the  laser  intensity  must  be  calibrated  in-situ  with  a  known  mixture, 
which  is  not  possible  in  many  engineering  applications  due  to  window  fouling,  vibrating  optics  and  beam  steering.  It 
has  been  demonstrated  that  dividing  the  second  harmonic  signal  by  the  first  harmonic  signal  removes  the  dependence 
on  laser  intensity,  thus  eliminating  the  need  for  intensity  calibrations  [1].  The  resulting  signal  is  referred  to  hereafter 
as  the  normalized  2f  signal.  A  lock-in  amplifier  is  used  to  isolate  the  harmonics,  filtering  noise  outside  of  the 
bandwidth  of  the  lock-in  amplifier,  and  shifting  the  detection  to  higher  frequencies  where  laser  and  system  noise  is 
decreased.  By  comparing  the  shape  and  magnitude  of  the  normalized  2f  signal  from  the  lock-in  detector  to  simulated 
signals,  thermodynamic  properties  such  as  pressure,  temperature  and  species  concentration  can  be  measured. 

II.  Experimental  Setup 

The  WMS  measurements  are  performed  at  the  Air  Force  Research  Laboratory  at  Edwards  Air  Force  Base  in 
California.  The  converging-diverging  nozzle  and  heating  system  was  mounted  inside  of  chamber  4  within  the 
Spacecraft  Propulsion  Laboratory.  The  laser  light  and  translation  stage  control  system  signals  were  passed  into 
chamber  4  via  fiber  optic  and  electrical  vacuum  feedthroughs  on  the  chamber. 

The  exhaust  plume  examined  in  this  study  is  composed  of  a  mixture  of  anhydrous  ammonia  and  water.  The  water 
concentration  of  the  plume  is  held  at  approximately  18%  water  by  mass.  The  gas  flow  is  preheated  by  two  circulation 
heaters  connected  in  series.  Liquid  water  from  a  syringe  pump,  is  injected  between  the  two  heating  stages,  exposing 
it  to  the  preheated  ammonia  flow.  The  syringe  pump  allows  the  water  flow  rate  to  be  precisely  controlled  and  varied. 
The  temperature  of  the  tube  where  the  water  enters  the  ammonia  flow  stream  is  monitored  by  a  type  K  thermocouple. 
To  adequately  vaporize  the  water,  water  is  not  injected  into  the  system  until  the  tube,  at  the  entry  point,  reaches  a 
minimum  temperature  of  100  °C.  To  minimize  system  exposure  to  the  corrosive  effects  of  ammonia,  the  system  is 
preheated  and  purged  with  nitrogen. 
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Immediately  following  the 
second  heating  stage,  the  gas 
mixture  flows  through  the 
converging -diverging  nozzle  seen 
in  Figure  1  to  create  the  plume. 
This  nozzle  configuration  results 
in  an  expansion  ratio  of  5.19  and 
an  exit  Mach  number  of  3.05.  By 
design,  this  nozzle  creates  a  highly 
under-expanded  plume  at  the 
operational  background  pressures 
of  550-720  mTorr. 


Three  tunable  diode  lasers  are 
used  to  probe  three  distinct 
absorption  features.  Laser  1 
(Eblana  EP1388-DM)  and  laser  2  (Eblana  EP1392-DM)  probe  H2O  absorption  features  at  1388  nm  and  1392  nm 
respectively.  Laser  3  (Toptica  LD- 1470-00 10-DFB-l)  probes  an  NH3  feature  at  1527nm.  The  transition  at  1388nm 
reaches  a  maximum  absorption  at  low  temperature,  while  the  transition  at  1392nm  reaches  a  maximum  absorption  at 
elevated  temperatures  (approximately  650K).  The  different  temperature  dependencies  of  the  H2O  absorption  features, 
results  in  an  absorbance  ratio  of  the  two  lines  that  is  sensitive  to  temperature,  but  also  results  in  low  peak  absorbance 
at  low  temperatures  for  the  feature  at  1391.7nm. 

A  National  Instruments  PXIe-6368  data  acquisition  card  provides  the  voltage  signal  used  to  modulate  the  laser 
injection  currents  for  WMS.  The  light  from  the  three  fiber-coupled  lasers  is  multiplexed  onto  a  single-mode  fiber  and 
passed  into  the  vacuum  chamber  through  an  optical  fiber  feed-through,  to  a  launch  optic  (Thorlabs  CFC-5X-C).  The 
beam  traverses  the  plume  perpendicular  to  the  flow  direction,  as  shown  in  Figure  1 .  The  beam  is  collected  by  a 
catch  collimator  (Oz  Optics  HPUCO-25-1300/1550-M-25AC-SP)  to  a  400-micron  diameter  multimode  fiber 
(Thorlabs  M45L01)  and  passes  back  through  the  optical  feed-through  to  an  InGaAs  photodetector  (Thorlabs 
PDA10CS).  The  signal  from  the  photodetector  is  monitored  by  the  PXIe-6368  data  acquisition  card  at  a  2MHz 
sampling  rate. 

The  launch  and  collection  optics  are  mounted  on  a  2 -axis  stage  that  can  be  translated  in  the  vertical  and  horizontal 
directions  through  the  use  of  two  stepper  motors.  The  translation  system  enables  measurements  from  -0.4  cm  to  +0.4 
cm  relative  to  the  plume  center  at  vertical  increments  of  0.1  mm,  resulting  in  81  total  measurement  locations  at  1  mm 
downstream  of  the  nozzle  exit. 


Gaseous 
NH  “ 


Liquid  H20 

I 


Heater 


Heater  2 


•  Measurement  Location 
%  Pitch  Optic 
■  Translation  Stage 


't? 


Figure  1:  Experimental  setup  of  NH3  and  H2O  injection  with  translation 
stage 


III.  Data  Processing 


A.  Residual  Amplitude  Modulation 

The  magnitude  of  the  absorption  signal  is  linearly  correlated  with  the  path  length  and  approximately  linearly 
correlated  with  pressure,  as  seen  in  Equation  1.  The  small  exit  diameter  of  the  nozzle  (4.57mm)  and  the  low  operating 
pressures  within  the  vacuum  chamber  (approximately  500  mTorr)  resulted  in  low  absorption  signals. 


Due  to  the  low  pressure  and  short  pathlength,  the  magnitude  of  the  absorption  signal  is  on  the  same  order  of 
magnitude  as  the  residual  amplitude  modulation  (RAM)  of  the  laser  [8],  [9].  RAM  is  caused  by  non-linear  intensity 
tuning  of  the  laser  system  that  creates  harmonic  signals  even  in  the  absence  of  absorption  (WMS  background).  It  is 
present  in  varying  degrees  in  all  tuned  laser  systems.  Figure  2  shows  the  normalized  2f  signal  of  the  transmitted 
intensity  (i.e.  after  absorption)  in  green  and  the  normalized  2f  signal  of  the  incident  intensity  (i.e.  without  absorption), 
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in  red.  Theoretically,  if  the  laser  tuning  is  exactly 
linear,  the  normalized  2f  of  the  incident  intensity 
would  be  constant  at  zero. 

In  Figure  2a  it  is  clear  that  the  non-linear  tuning 
creates  a  strong  background  in  the  normalized  2f 
signal  of  the  transmitted  intensity.  The  incident 
intensity  is  recorded  concurrently  with  the 
transmitted  intensity;  this  allows  the  second 
harmonic  components  of  the  non-absorbing 
incident  intensity  signal  to  be  subtracted  from  the 
second  harmonic  components  of  the  absorbing 
transmitted  intensity  signal  according  to  Equation 
2.  The  parameters  X2f,  Y2f  are  the  X  and  Y 
components  of  the  signals  from  the  digital  lockin  at 
2f.  Rif  refers  to  the  root  mean  square  of  the  X  and 
Y  parameters  at  If.  Finally,  the  components 
denoted  ‘raw’  are  from  the  transmitted  intensity 
signal  before  RAM  subtraction,  while  the 
components  denoted  ‘bg’  or  background  are  from 
the  incident  intensity  signal.  With  this  procedure, 
we  can  isolate  the  components  of  the  transmitted 
intensity  signal  caused  by  absorption  The 
background- subtracted  signal  is  shown  in  Figure 
2b.  This  method  is  applied  to  all  measured  signals 
for  all  three  lasers.  The  1392nm  feature,  shown  in 
Figure  2,  is  the  weakest  absorption  signal  (9E-4 
absorbance,  or  0.09%  absorption),  and  therefore  the 
most  sensitive  to  the  background  signal  and 
subtraction. 


Figure  2:  Figure  2a:  Normalized  2f  signal  of  transmitted 
intensity  with  absorption  before  RAM  subtraction  (green), 
Normalized  2f  signal  of  incident  intensity  without  absorption, 
showing  the  RAM  contribution  (red).  Figure  2b:  Normalized 
2f  signal  with  RAM  subtraction. 


B.  Averaging 

For  these  measurements,  we  perform  scanned- wavelength  WMS,  where  the  high  frequency  wavelength 
modulation  is  combined  with  an  additional  slow  sine  wave  modulation  (100Hz)  to  scan  the  average  wavelength  of  the 
laser  across  the  entire  absorption  feature  [10].  The  basic  time  resolution  of  the  measurement  is  therefore  200  Hz,  as 
we  scan  the  full  normalized  2f  signal  for  each  feature  twice  per  cycle  of  the  slow  modulation. 

Figure  3  shows  a  single  sweep,  and  a  0.5s  average  of  the  normalized  2f  signals.  Though  the  single  and  0.5s 
averages  are  similar  for  the  stronger  features  probed  by  lasers  1  and  3,  the  0.5s  average  significantly  improves  data 
analysis  for  the  weak  1392nm  feature.  For  the  results  reported  in  the  remainder  of  the  paper,  0.5s  averages  are  used. 


Figure  3:  Raw  data  (red)  versus  time  averaged  data  (blue)  for  each  absorption  feature.  Figure  3a:  1388nm  H2O  feature, 
Figure  3b:  1392nm  H2O  feature,  Figure  3c:  NH3  feature. 
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C.  Absorbance  Area 

We  utilize  a  fitting  procedure  similar  to  that  described 
by  [11]  to  calculate  an  integrated  absorbance  area  for  each 
absorption  feature  [6],  [3],  [5].  Using  the  integrated 
absorbance  area,  rather  than  the  peak  normalized  2f  signal, 
simplifies  the  interpretation  of  the  measured  data  by 
removing  the  dependence  of  the  signal  on  line  broadening. 
An  example  of  a  fit  produced  with  this  procedure  is  shown 
in  Figure  4.  The  good  agreement  between  the  fit  and  data 
ensures  that  the  absorbance  area  calculated  by  the  fit  is 
representative  of  the  absorbance  area  of  the  measured 
absorption  feature. 


Figure  4:  Fitted  normalized  signals  (black)  versus 
The  fitting  procedure  is  repeated  for  each  measurement  ^me  averaged  normalized  signals  (blue)  for  1527nm 
location  within  the  81 -point  vertical  scan  at  1  mm  NF[3  feature  at  1mm  downstream  of  nozzle  exit, 
downstream  of  the  nozzle  exit.  The  resulting  areas  for  all 

three  absorption  features  at  each  point  within  the  scan  are  displayed  in  Figure  5  as  a  function  of  vertical  position  of 
the  beam  within  the  plume.  The  absorbance  areas  do  not  decrease  to  zero  at  the  edges  of  the  plume  due  to  residual 
water  vapor  and  ammonia  within  the  vacuum  chamber.  To  account  for  the  residual  signal,  absorbance  areas  for  all 
three  absorption  features  are  measured  outside  of  the  plume  and  are  subtracted  to  isolate  absorbance  within  the  plume 
for  temperature  and  species  concentration  calculations. 


1388  nm  1392nm  1527nm 


Figure  5:  Absorbance  areas  as  a  function  of  vertical  position  in  plume,  a)  Absorbance  area  of  1388nm  FFO  feature, 
b)  Absorbance  area  of  1392nm  H2O  feature,  c)  Absorbance  area  of  1527nm  NH3  feature. 

IV.  Path  Averaged  Results 


A.  Temperature 

Two-line  thermometry,  described  by  [12],  is  implemented  to  calculate  the  path  averaged  temperature  within  the 
plume.  By  taking  the  ratio  of  the  calculated  absorbance  areas  from  the  two  H2O  absorption  features,  the  measurement 
dependence  on  pressure,  path  length  and  species  concentration  is  removed,  as  seen  in  Equation  3.  Under  these 
circumstances,  the  ratio  of  the  absorbance  areas,  A1388  and  A1392,  simplifies  to  the  ratio  of  the  linestrengths  of  the  H2O 
absorption  features,  Si38s(T)  and  Si392(T),  which  are  solely  a  function  of  temperature.  Two-line  thermometry  is 
especially  useful  for  plume  measurements,  as  the  plume  width  (L)  and  pressure  distributions  are  not  precisely  known 
at  each  measurement  location. 

^1388  _  ^1388 (TV^H2Q^Ptotal  _  $1388  CO  (3) 

^1392  Sl392(T)XH2oLPtotal  *^1 3  92  CO 
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Using  this  approach,  the  path-averaged  temperature  at  each  point  within  the  plume  at  1mm  downstream  of  the 
nozzle  exit  plane  is  calculated  and  shown  in  Figure  6.  At  the  edges  of  the  plume  the  absorbance  signal  from  residual 
H2O  and  NH3  make  up  a  larger  percentage  of  the  total  measured  signal,  as  the  absorbance  signal  contributed  from  the 
plume  decreases.  This  causes  the  uncertainty  associated  with  the  calculated  background  absorbance  to  be  magnified 
at  the  edges  of  the  plume  as  the  background  becomes  a  larger  percentage  of  the  signal. 


Figure  6:  Path  averaged  temperature  profile  as  a  function  of  vertical 
position  within  the  plume. 


B.  Species  Concentration 

Common  methods  for  calculating  species  concentration  in  the  plume  are  not  effective  for  this  application  as  the 
pathlength  of  the  laser  through  the  plume  and  the  pressure  distribution  are  unknown.  Here  we  introduce  a  new  method 
of  species  concentration  measurement  that  is  possible  when  one  is  measuring  all  of  the  species  present  in  the  flow.  To 
begin,  we  calculate  the  ratio  of  the  absorbance  areas  of  one  of  the  water  absorbance  features  and  the  ammonia 
absorbance  feature.  Both  absorbance  areas  were  recorded  at  the  same  path  length  (L)  and  pressure  (P total),  therefore 
the  dependence  on  these  parameters  is  removed  (Equation  4). 

In  this  system  we  assume  there  are  no  reactions  taking  place,  therefore  the  only  molecules  within  the  system  are 
the  ammonia  and  water  vapor.  Under  this  situation,  the  sum  of  the  mole  fractions  of  ammonia  and  water  are  equal  to 
one. 


^1388  $1388(T)XH20LPt 


total 


>1388 


cm. 


H20 


^1527  $ 1527  (T)XNH3LPtotai  $1527(T)XnH3 

N 

Xi  =  1  ->  XNH3  +  XH2  0  ~  1 


i- 


XnH3  — 


J1388 


(T)A 


1527 


^1388*^1527  CO  +  ^1527*^1388  CO 


(4) 

(5) 

(6) 
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By  combining  Equations  4  and  5  and  solving  for  the  mole  fraction  of  ammonia  we  arrive  at  Equation  6.  Since  the 
path  averaged  temperature  has  already  been  calculated,  the  right  side  of  this  equation  is  entirely  known,  enabling  the 
calculation  of  the  path-averaged  mole  fractions  of  the  absorbing  species  at  each  location  in  the  scan  (see  Figure  7). 
Since  the  vacuum  chamber  was  evacuated  prior  to  trials,  we  expect  the  mole  fraction  of  NH3  and  H2O  to  be 
approximately  constant  throughout  the  plume.  The  variation  in  mole  fractions  in  Figure  7  is  likely  due  to  incomplete 
mixing  of  the  NH3  and  H2O  prior  to  the  nozzle,  or  non-linear  temperature  effects  resulting  from  path  averaged 
measurements  over  a  range  of  temperatures.  The  large  uncertainties  on  the  edges  of  the  plume  is  due  to  temperature 
uncertainty,  which  is  reflected  in  line  strength  uncertainty  in  Equation  6. 
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Figure  7:  Water  and  ammonia  mole  fraction  distributions  in  plume.  Dashed  lines  represent  the  mean  mole 
fraction  of  the  data.  Solid  black  lines  represent  the  expected  mole  fraction  based  on  the  flow  rates  of  H2O  and 
NH3  into  the  chamber. 


The  dashed  lines  in  Figure  7  represent  the  average  mole  fractions  of  H2O  and  NH3  from  the  data,  0.167+0.013  and 
0.833+0.013  respectively.  The  solid  black  lines  represent  the  expected  mole  fraction  based  on  the  flow  rates  of  H2O 
and  NH3  into  the  chamber,  0.172  and  0.828  respectively.  The  expected  mole  fractions  fall  well  within  one  standard 
deviation  of  the  calculated  values.  The  agreement  between  expected  values  and  measured  values  suggests  that  the 
proposed  species  concentration  technique  is  valid  for  this  application  and  that  all  significant  constituents  of  the  flow 
are  accounted  for. 

Discrepancies  in  the  measured  mole  fractions  versus  the  expected  mole  fractions  could  be  due  to  reactions  taking 
place  within  the  system,  which  we  did  not  account  for  by  assuming  the  only  constituents  in  the  flow  were  H2O  and 
NH3,  or  from  water  vapor  condensing  out  of  the  flow  due  to  rapid  temperature  drops  from  the  rapid  expansion  into  the 
vacuum  chamber.  At  the  elevated  temperatures,  such  as  the  temperature  within  the  heaters  (approximately  800K), 
NH3  will  partially  to  dissociate  into  N2  and  H2  as  in  Equation  7.  While  the  resident  time  the  ammonia  was  kept  at 
elevated  temperatures  was  short,  this  reaction  likely  decreased  the  mole  fraction  of  the  NH3  present  in  the  system. 
This  would  also  decrease  the  H2O  mole  fraction  since  the  total  number  of  moles  in  the  system  is  increased  by  this 
reaction. 


2NH3  <=>  N2  +  3 H2 


(7) 


Although  condensation  is  not  observed  within  the  system  during  testing,  it  is  possible  that  the  H2O  vapor  is 
condensing  in  the  flow,  particularly  at  the  edges  of  the  plume  where  the  temperatures  are  the  lowest.  Condensation 
would  cause  the  vapor  mole  fraction  of  H2O  to  decrease,  while  increasing  the  vapor  mole  fraction  of  NH3.  This  is 
consistent  with  the  calculated  mole  fractions  being  lower  than  the  expected  values  in  Figure  7. 
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V.  Abel  Transform  Results 


Given  the  radially  symmetric  nozzle  geometry  and  the  symmetry  of  the  path  averaged  measurement  results  about 
the  center  axis  of  the  nozzle  exit,  we  assume  that  the  plume  is  approximately  axisymmetric.  This  allows  the  path 
averaged  integrated  absorbance  measurements  (. Ax(y ))  from  Figure  5  to  be  converted  into  integrated  absorbance  as  a 
function  of  radius  (Ax(r)),  by  the  use  of  an  inverse  Abel  transform  (Equation  8).  This  method  is  described  in  depth  in 


[13]. 


A\(r) 


1  rR  dAx  (y)  dy 
n)y  dy  Jy2  _  r2 


(8) 


We  can  reconstruct  the  radial  temperature 
profile  within  the  plume  using  two -line 
thermometry  and  the  radial  absorbance  area 
distributions  of  the  two  H2O  features  (Figure  8). 
From  isentropic  flow  calculations  for  pure  NH3, 
we  estimate  that  the  temperature  at  the  core  of  the 
plume  at  the  exit  of  the  nozzle  should  be 
approximately  350K.  As  our  measurement  is 
downstream  of  the  exit  plane,  it  is  logical  that  the 
measured  temperatures  are  lower  than  this 
predicted  exit  temperature.  Computational  fluid 
dynamics  models  would  be  required  to  determine 
the  expected  temperatures  at  the  measurement 
location.  The  high  temperature  core  and  the  low 
temperature  edges  of  the  plume  are  consistent 
with  the  expected  nearfield  under-expanded 
structure  [14],  [15].  As  the  optics  are  translated 
past  to  the  plume  boundary,  we  see  the 
temperature  return  to  near  atmospheric 
temperature. 


Figure  8:  Radial  temperature  distribution  from  inverse  Abel 
transform 


VI.  Conclusion 

We  resolve  path-averaged  and  radial  profiles  of  temperature  within  the  plume  of  a  small  diameter  (<5mm) 
converging-diverging  nozzle  expanding  into  a  vacuum  chamber  at  Edwards  Air  Force  Base.  A  new  method  of 
calculating  species  mole  fraction  was  proposed,  for  systems  in  which  pressure  and  path  length  are  uncertain,  and  for 
which  all  species  are  being  measured.  While  the  species  mole  fraction  measurement  technique  includes  some  obvious 
sources  of  uncertainty,  the  uncertainties  associated  with  estimating  pressure  and  path  length  in  a  system  such  as  this 
would  be  far  greater.  The  calculated  mole  fractions  compare  well  with  expected  values,  indicating  this  mole  fraction 
calculation  technique  is  valid.  Temperature  results  from  path  averaged  and  radially  distributions  agree  with  theoretical 
predictions,  but  a  computational  fluid  dynamics  model  would  be  needed  to  fully  validate  this.  The  positive  results 
from  this  study  indicate  that  WMS  could  be  applied  to  operational  propulsion  systems  to  gain  temporally  and  spatially 
resolved  measurements  within  the  thruster  plumes. 
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Wavelength  modulation  spectroscopy 
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Absorption  Spectroscopy 


■  Molecules  absorb  light  at  discrete  wavelengths 
resonant  with  molecular  transitions 

■  The  amount  of  light  absorbed  is  related  to: 

a  (A)  =  S(T)PXL$> 

-  Line-strength:  s(T) 

-  Pressure:  P 

-  Mole  fraction:  X 

-  Pathlength:  L 

-  Line-shape  function:  4 

■  Parameters  S(T)  and  <P  are  unique  to  each 
transition 

-  HITRAN2012[3] 

-  Experimental  validations  [4,5] 


University  of  Colorado 
Boulder 


Distribution  Statement  A:  Approved  for  Public  Release;  Distribution  is  Unlimited.  PA#  17443 


Tunable  Diode  Laser  Absorption  Spectroscopy 

(TDLAS) 


Ramp 

/VI  IS 

Laser 


Gas 

Sample 


->L 

Detector 


University  of  Colorado 
Boulder 


Distribution  Statement  A:  Approved  for  Public  Release;  Distribution  is  Unlimited.  PA#  17443 


Tunable  Diode  Laser  Absorption  Spectroscopy 

(TDLAS) 


Ramp 

/VI  IS 

Laser 


Gas 

Sample 


->L 

Detector 


q(A)  =  S{T)PXL<5>  =  -ln(j-)  < 


Wavelength 


Wavelength  Modulation  Spectroscopy 


Ramp 


Sine  of  frequency  f 


University  of  Colorado 
Boulder 


Distribution  Statement  A:  Approved  for  Public  Release;  Distribution  is  Unlimited.  PA#  17443 


9 


Wavelength  Modulation  Spectroscopy 


Ramp 


Gas 

Sample 


Detector 


Sine  of  frequency  f 


Time 


University  of  Colorado 
Boulder 


Distribution  Statement  A:  Approved  for  Public  Release;  Distribution  is  Unlimited.  PA#  17443 


Wavelength  Modulation  Spectroscopy 


Ramp 


Sine  of  frequency  f 


1000- 


*• 


Time 


University  of  Colorado 
Boulder 


Distribution  Statement  A:  Approved  for  Public  Release;  Distribution  is  Unlimited.  PA#  17443 


11 


Magnitude 


Wavelength  Modulation  Spectroscopy 


Ramp 


Gas 

Sample 


Detector 


Sine  of  frequency  f 


Fast  Fourier  Transform 
1/  , _ 2/ 


Frequency  (kHz) 


1000- 


Time 


University  of  Colorado 
Boulder 


Distribution  Statement  A:  Approved  for  Public  Release;  Distribution  is  Unlimited.  PA#  17443 


12 


Magnitude 


Wavelength  Modulation  Spectroscopy 


1031 


Frequency  (kHz) 


[6] 


University  of  Colorado 
Boulder 


Distribution  Statement  A:  Approved  for  Public  Release;  Distribution  is  Unlimited.  PA#  17443 


3 


Wavelength  Modulation  Spectroscopy 


University  of  Colorado 
Boulder 


Distribution  Statement  A:  Approved  for  Public  Release;  Distribution  is  Unlimited.  PA#  17443 


14 


2f  Mag.  If  Mag. 


Wavelength  Modulation  Spectroscopy 


Normalizing  the  2f  signal  by 
the  If  signal  removes 
dependence  on  laser  intensity 

Important  in  harsh 
environments 

-  Window  fouling 

-  Vibrating  optics 


Vibrating  Optics 


CD 

73 

3 

'c 

S> 

05 


ft 


0.00- 


CD 

73 

3 

'c 

s> 

co 


CD 

TD 


co 


ft 


0.0- 


0.00- 


0.0  0.1  0.2 


0.3  0.4  0.5 

Time  (s) 


j  Y 

2f  Magnitude 

Y 

Y . 

n  1  i 

i  |  i  |  i 

1  1  1  1  1 

1  1  1  1 

:  ir 

-  If  Magnitude 

V 

Y 

n  1  i 

i  |  i  |  i 

1  1  1  1  1 

1  1  1  1 

:  2f/1 

- 1 - 1 

f 

- 1 - 1 - 1 - 1 - 1 

- 1 - 1 - 1 - 1 - 1 

1 - 1 - 1 - 1 

[6,7] 


0.6  0.7 


University  of  Colorado 
Boulder 


Distribution  Statement  A:  Approved  for  Public  Release;  Distribution  is  Unlimited.  PA#  17443 


15 


Data  Processing 


Distribution  Statement  A:  Approved  for 
Public  Release;  Distribution  is  Unlimited. 
PA#  17443 


Wavelength 

i  i  i  1 1 1 

V\/ .  i; 

\  Absorbance 

J  V  vs  Time  a(t) 

- 1 - 1 - 1 - 1 - 1 - “ - 1  . 

-|-jme  Wavelength 

Incident  Intensity 

i 

&  l0(t)e 

Q 

II 

ansmitted  Intensity 

i  i  i  i  i 

i  i  i  i  i  i 

Time 

p  i  i  i  i  i  i 

Time 

it 

Run  resulting  l(t)  «  -E2SI—  A 

through  lockin  to  generate  f  -  O  A 

a  simulated  signal.  Compare  |  -  /\J/  W\ 

[8,9] 


Temperature  Measurement 


Temperature  can  be  calculated  from  the  ratio  of  the 
absorbance  areas  of  two  well  chosen  absorption  features 


Absorbance  area  ratio  =  S,  }P XL 

S2(T)PXL 
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Absorbance 


Temperature  Measurement 
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System  Design 


■  Testing  performed  at  the  AFRL  at 
Edwards  Air  Force  Base 

-  Vacuum  chamber  4  within  the 
Spacecraft  Propulsion  Laboratory 

■  System  developed  by  the  AFRL 

-  Exhaust  products:  NH3  and  H20 

-  Expansion  ratio:  5.19 

-  Exit  Mach  number:  3.1 

■  LabView  controlled  translation 
stage  mounted  downstream 
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Spectroscopy  System 


■  Spectroscopy  system 
developed  at  CU  Boulder 

-  Verification  of  NH3and  H20  line 
parameters  in  furnace 

-  Characterized  laser  operating 
conditions 

-  Developed  simulation 
procedures 

■  Measurements 

-  Temperature  distributions 

-  H20  concentration 

-  NH3  concentration 

■  Three  Near-IR  lasers 

-  1388nm  (H20) 

-  1392nm  (H20) 

-  1527nm  (NH3) 
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Experimental  Setup 


Vertical  scan  Digital  lock-in 
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Scanning  Path 


x  =  1  mm  (from  exit  plane  of  nozzle),  81  point  vertical  scan,  Ay  =  0.1  mm 
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Scanning  Path 


x  =  1  mm  (from  exit  plane  of  nozzle),  81  point  vertical  scan,  Ay  =  0.1  mm 
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Data  Processing 


Water  absorption  features 


Typical  fit  for  1388nm  signal  (H20) 


Typical  fit  for  1392nm  signal  (H20) 


Good  agreement  between 
the  simulated  signal  and 
the  fit  indicate  correct 
absorbance  parameters 


[8,9] 
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Absorbance  Area 


H20  feature  at 
1388nm 


Data  recorded  at  1  mm  from  exit  plane  of  thruster 
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Absorbance  Area 


H20  feature  at 
1388nm 


H20  feature  at 
1392nm 


Absorbance  Area  (THz)  Absorbance  Area  (THz) 


Data  recorded  at  1  mm  from  exit  plane  of  thruster 
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Absorbance  Area 


H20  feature  at 
1388nm 


H20  feature  at 
1392nm 


IMH3  feature  at 
1527nm 


Absorbance  Area  (THz)  Absorbance  Area  (THz)  Absorbance  Area  (THz) 


Data  recorded  at  1  mm  from  exit  plane  of  thruster 
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Vertical  Position  (cm) 


Absorbance 


Temperature  Measurement 


Temperature  can  be  calculated  from  the  ratio  of  the 
areas  for  1388nm  and  1392nm  absorption  features 


Absorbance  area  ratio  = 


S(1388y*jft 

5(1392)J»rtl 


5(1388) 

5(1392) 
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Temperature  Measurement 


270  280  290  300  310  320 

Temperature  (K) 
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Vertical  Position  (cm) 


Inverse  Abel  Transform 


■  Apply  inverse  Abel  transform  to  path 
averaged  measurements 

-  Assuming  axisymmetric  flow 

-  Calculate  radial  distributions 

-  Use  “Fourier  Method”  to  approximate 
F(y)  as  a  cosine  expansion 


fRdF(y)_dy__ 
y  jy2  ~  r2 


[11,12] 
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Inverse  Abel  Transform 


Absorbance  Area  of  1388  Feature 


[11,12] 


University  of  Colorado 
Boulder 


Distribution  Statement  A:  Approved  for  Public  Release;  Distribution  is  Unlimited.  PA#  17443 


36 


Inverse  Abel  Transform 


Absorbance  Area  of  1388  Feature 


y,  r 


y,  r  (cm) 


[11,12] 
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Inverse  Abel  Transform 


Absorbance  Area  of  1388  Feature 


y,  r 


y,  r  (cm) 


[11,12] 
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Inverse  Abel  Transform 


■  Isentropic  calculations  predict  core 
temperature  of  353K  at  outlet 

-  Results  are  1  mm  downstream  of  outlet 

-  Magnitude  of  temperature  drop  agrees 
with  Naik  et  al.  (2009)  and  Woodmansee 
et  al.  (2004) 

■  Transverse  temperature  trends  are 
consistent  with  Woodmansee  et  al. 
(2004) 

-  Exit  plane  temperature  measurements 

■  Lower  temperature  regions  coincide 
with  edges  of  nozzle 


[11,12] 
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Species  Concentration 


■  Traditionally  to  calculate  concentration  (X)  from  absorbance  area  we  must 
know  pressure  and  path  length 

-  Pressure  distribution  is  unknown  at  measurement  location 

-  Edges  of  the  plume  are  not  distinctly  defined 


■  Computational  models  could  be  used 

-  Time  intensive 

-  Could  introduce  large  errors  if  not  correct 

■  Simulating  pressure  distributions  adds  complexity  to  absorbance  area 
simulation 
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Species  Concentration 


New  technique  similar  to  the  temperature  measurement  technique  shown 
earlier 

-  Take  the  ratio  of  absorbance  areas  of  the  NH3  feature  and  the  1388  nm  H20  feature 


^1388  _  ‘-*1388  (TV^H2Q^^ebL  _  ‘-*13  88  (T)XH20 
^1527  ‘$1527  (T)XNH3 

^1527(^)^NH3 

Assume  NH3  and  H20  are  the  only  significant  species 

i  N 

'  =  1  XNH3  +  XH2 o  —  1 


ZiV 

*  = 

i=l 


(1) 


(2) 


-  Solve  for  XNH3  by  combining  Equation  1  and  2 


Xnh3  ~ 


•$1388  00^152  7 


^1388*^1527  (T)  +  ^1527*^1388  (T) 
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Species  Concentration 


-  Mean  mole  fractions  from  data 

-  XH20  =  0.167±0.013,  XNH3  =  0.833±0.013  (dashed  lines) 

■  Expected  mole  fractions  based  on  measured  flow  rates 

-  XH20  =  0.172,  XNH3  =  0.828  (solid  lines) 

■  Uncertainty  on  edges  due  to  temperature  uncertainty  in  S(T) 
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WMS  as  a  Diagnostic 


Benefits: 

-  Small  robust  system 

-  Relatively  inexpensive 

■  Easy  alignment 

■  Calibration  free 

-  Measurements  of  temperature  and  mole 
fraction 

■  Low  power 

-  Fast  acquisition  times 

-  Non-intrusive 

Drawbacks: 

■  Path  averaged  1 D  measurements 

■  Velocity  measurements  are  challenging 
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Planar  Laser  Induced 
Fluorescence  (PLIF) 


WMS  as  a  Diagnostic 


Benefits: 

-  Small  robust  system 

-  Relatively  inexpensive 

■  Easy  alignment 

■  Calibration  free 

-  Measurements  of  temperature  and  mole 
fraction 

■  Low  power 

-  Fast  acquisition  times 

-  Non-intrusive 

Drawbacks: 

■  Path  averaged  1 D  measurements 

■  Velocity  measurements  are  challenging 


Wavelength  Modulation 
Spectroscopy  (WMS) 
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Conclusions  and  Acknowledgements 


Conclusions: 

■  Spatial  distributions  of  temperature  were  calculated  within  the  plume 

■  Radial  temperature  distributions  were  calculated  using  Abel  inversions 

-  Temperature  distributions  consistent  with  previous  experimental  results  [13,14] 

■  New  method  for  calculating  species  concentration  was  proposed 

-  Concentration  measurements  compare  well  with  expected  results 

■  Results  indicate  WMS  is  a  promising  technique  for  microthruster  diagnostics 
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Data  Processing 


Data  sets  were  averaged 
for  0.5  seconds  to 
increase  signal  to  noise 
ratio  (SNR) 
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Residual  Amplitude  Modulation  (RAM) 


■  Non-linearity  in  the  laser  2 
was  substantial 

■  l0  signal  was  used  for 
background  subtraction  to 
account  for  this 


Wavelength  Modulation  Spectroscopy 


X  Component: 
Y  Component: 

Magnitude: 


A  sin (fx  +  <j>)  cos (fx)  =  y  ( sin(0)  +  sin (2  fx  +  b)) 
A  sin  (fx  +  <p)  sin  (fx)  =  j  (  cos(^)  —  cos(2/x  +  b)) 

i /  (4  sin(|?i))2  +  (4  cos(^))2  =  a/2 
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Data  Processing 
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Inverse  Abel  Transform 


Radial  Position  (cm) 
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